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Abstract 



Signatures of neutrino and antineutrino signals from stellar collapse in 
heavy water detectors are contrasted with those in water detectors. The ef- 
fects of mixing, especially due to the highly dense matter in the supernova 
core, are studied. The mixing parameters used are those sets allowed by cur- 
rent understanding of available neutrino data: from solar, atmospheric and 
laboratory neutrino experiments. Signals at a heavy water detector, espe- 
cially the dominant charged current reactions on deuteron, are very sensitive 
to some of these sets of allowed mixing parameters. Theoretical uncertainties 
on supernova neutrino spectra notwithstanding, a combination of supernova 
measurements with water and heavy water detectors may be able to distin- 
guish many of these mixing possibilities and thus help in ruling out many of 
them. 
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I. INTRODUCTION 



Neutrinos from stellar collapse have so far been detected only from the supernova 
SN1987a in the Large Magellanic Cloud Jl|. The inital observations of the neutrino and 
antineutrino events by the Kamiokande @ and 1MB || collaborations were the subject of 
detailed analysis by several authors fl-|T0"f immediately following the event. The analyses 
confirmed the qualitative features of core collapse and subsequent neutrino emission. 

The effect of non-zero neutrino masses and mixing on supernova signals was first analysed 
in general by Kuo and Pantaleone [|llj] . Recently, several authors have looked at the possible 
signatures of neutrinos and antineutrinos from supernova collapse in realistic scenarios. 
Dighe and Smirnov Jl2| have looked at the problem of reconstruction of the neutrino mass 
spectrum in a three-flavour scenario. These authors as also Chiu and Kuo |TB| have compared 
the signatures in the standard mass hierarchy and inverted mass hierarchy. While these 
papers incorporate the constraints from solar and atmospheric neutrino observations, the 
important question of the mass limits that may be obtained from the observation of time 



delay has been analysed by Beacom and Vogel fll4]]l5l and by Choubey and Kar JT6[ (see 
also the review by Vogel [[H]]). For a recent review which also discusses aspects of locating 



a supernova by its neutrinos in advance of optical observation, see Ref. [|18 



In this paper we apply the analysis of neutrinos from stellar collapse presented in 
Refs. [0,0] to heavy water detectors. We had previously discussed in detail the signa- 



tures, in a water Cerenkov detector, of neutrinos and antineutrinos from stellar collapse in 
both 3- and 4-flavour mixing scenarios. The 4-flavour analysis was motivated by lsnd data 



2lfl . The analysis was confined to Type II supernovae (which occur when the initial mass 
of the star is larger than 8 solar masses). The choices of mixing parameters used were con- 
sistent with available data on solar, atmospheric and laboratory neutrino experiments. It 
turns out that different choices of (allowed) mixing parameters lead to drastically different 
supernova signals at water detectors. While certain features are common to both the 3- 
and 4-flavour analyses, there are important differences that may (depending on the mixing 
parameters) be able to distinguish the number of flavours. We will summarise the salient 
features of this analysis below. 

The dominant contribution is from the charged current (CC) scattering of V e on protons 
in water. Mixing can increase the number of high energy events in this channel. The 
most dramatic effect of 3-flavour mixing is to produce a sharp increase in the CC events 



involving oxygen targets [H|22] for most choices of mixing parameters. These will show 



up as a marked increase in the number of events in the backward direction with respect 
to the forward peaked events involving electrons as targets (more than 90% of which lie 
in a 10° forward cone with respect to the supernova direction for neutrinos with energies 
E v Z 8 MeV), both showing up over the mostly isotropic CC events on protons. In the 
absence of such mixing, there will be only a few events due to CC scattering on oxygen 
targets. These events involving oxygen targets will be seen in a heavy water detector as 
well; realistically however, there will be fewer such events due to the considerably smaller 
size of the heavy water detector at the Sudbury Neutrino Observatory (sno) as compared 
to the water detector at SuperKamiokande (SuperK). 

When 4-flavour mixing is considered, the analysis becomes obviously more complex; now, 
the increase due to oxygen events will be visible only for some of the allowed values of the 
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parameters. Furthermore, there is no set of allowed parameters in both the 3- and 4-flavour 
cases for which the signals in a water detector will be able to distinguish between adiabatic 
and non-adiabatic neutrino propagation in the core of the supernova. This is an important 
issue since it can place a lower bound on the (13) mixing angle, which is currently bounded 
by chooz at the upper end in the case of 3-flavour mixing. (The only non-zero value for this 
angle so far comes from the lsnd data which can only be analysed in a 4-flavour frame- work). 
A zero value for this mixing angle will decouple the two sets of oscillations, v e — > and 
Vn —>■ v T i and allow for a simple 2-flavour analysis separately of the solar and atmospheric 
neutrino data. An exactly zero value of this angle will also render irrelevant CP violating 
phases in the problem. 

There are also CC and neutral current (NC) events due to elastic scattering with electrons 
in water/heavy water. These interactions are the same in both detectors and have been 



discussed in detail in the earlier analysis on water detectors |20[ . 

In this paper we focus on the possible signatures of a supernova collapse in a heavy 
water detector. The question is of practical interest since SNO is operating for more than 
a year now. It turns out that a combination of measurements in water and heavy water 
detectors has much better discriminatory power than either of them individually. Hence such 
observations of supernova neutrinos may be a good signal to rule out some of the currently 
allowed parameter space in the neutrino mixing angles. (These signals are however not very 
sensitive to neutrino mass squared differences.) 

The most interesting signals to study at a heavy water detector are the events from CC 
interactions of both v e and V e on deuterons. This is the most dominant signal in contrast to 
the dominant CC interaction of V e alone on free protons in a water detector (all are typically 
about two orders of magnitude larger than those from oxygen or electron interactions). We 
will focus mostly on these events in this paper, besides making a few remarks on the NC 
events on deuteron that can also easily be measured at a heavy water detector. These are 
therefore the "new" signals in a heavy water detector that would not be observable in (or 
will be different from) a water detector. 

As before the analysis is done assuming the standard mass hierarchy necessitated by the 
solar and atmospheric neutrino observations |23| . We also impose all the known constraints 
on the mixing parameters and mass-squared differences including those constraints from 
laboratory experiments pT[ , p4]| . The purpose of the calculation is to see whether different 
mixing scenarios give significantly different signals in the detector. 

In Sect. II we briefly outline the frame- work including the mixing matrix as also matter 
effects on mixing. We then use this to obtain expressions for the u e , V e and v^ T , V^ T fluxes 
reaching the detector in both 3- and 4-flavour mixing schemes. While discussions may be 



found in Ref. [19, 20 1 we reproduce the relevant details to keep this paper self-contained. 



In Sect. Ill we list the different interaction processes relevant to both water and heavy 
water detectors, in particular, those with electrons (positrons) in the final state, as well 
as NC events on deuteron. The details of the supernova model used in the calculation as 
well as the allowed values of neutrino mixing parameters (from already existing data) using 
which the effects of mixing are computed are also listed here. Results showing the effects of 
mixing both for the spectrum and the integrated number of events are presented in Sect. IV. 
Sensitivity to the supernova model parameters is discussed. The NC events on deuterons 
in a heavy water detector are also briefly discussed here. All these results are for a (2 + 2) 
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mass hierarchy scheme in the case of 4 flavour mixing. In Sect. V, we briefly discuss the 
other possible scenario, that of the (3 + 1) scheme in 4-flavours. We present a summary and 
discussions in Sect. VI. 



II. NEUTRINO MIXING AND MATTER EFFECTS 

We briefly review mixing among three and four flavours of neutrinos (or antineutrinos) 
and compute the neutrino (antineutrino) survival and conversion probabilities. These are 
given in more detail in Refs. [|19|,[2(J. The supernova neutrinos are produced mostly in the 
core, where the matter density is very high. Hence matter effects on the propagation are 
important. 

The hierarchy of mass eigenstates is shown in Fig. |]. The 3-flavour case is shown in 
Fig. [TJa; there are essentially two scales corresponding to the solution of solar and atmospheric 
neutrino problems assuming neutrino oscillations. In the case of four flavours (necessitated 
by the non-zero result of lsnd ||21||), the analysis is more complicated since there is an 
additional sterile neutrino. One possible mass hierarchy here is a (2 + 2) scenario [25] as 
shown in Fig. [l|b. The lsnd result implies the existence of a mass scale in the range of 
0.1 eV 2 to 1 eV 2 . We choose two doublets separated by this mass scale. The intra-doublet 
separation in the lower doublet corresponds to the solar neutrino scale ~ 10~ 5 eV 2 and 
that in the upper one to the atmospheric neutrino scale ~ 10 -3 eV 2 . Yet another possibility 



is the so-called (3 + 1) scheme |26|j27|1 , where the 3-flavour scheme is extended by adding 



a heavy sterile neutrino as the heaviest mass eigenstate, with a mass separation from the 
other three states as required by lsnd. We discuss this scheme separately later. In what 
follows, therefore, 4-flavour mixing always refers to the (2 + 2) scheme. 

The mixing matrix, U, which relates the flavour and mass eigenstates in the above 
scenarios has three angles in the 3-flavour case and six angles in the 4-flavour one. We shall 
ignore the CP violating phases. Then, the mixing matrix can be parametrised in the case 
of 3-flavours as, 



We v r\ = U X [V X V 2 U 3 ] , (1) 

where T stands for transpose. Here U is parametrised by considering rotations of mass 
eigenstates taken two at a time: 



U = U 23 x U 13 x U 12 



= x [7 x U w , (2) 

where Uij corresponds to the rotation of mass eigenstates |z) and As is the convention, 
we denote the mixing angle relevant to the solar neutrino problem by u and that relevant to 
the atmospheric neutrino problem by ip in the case of 3 flavours. The (13) angle is then 
small, due to the chooz result which translates to the constraint on the (13) mixing 
angle, sin0 = e < e , where e = 0.16. In the case of 4-flavours, we choose to work in the 
(2 + 2) scheme: 

[u e u s u T ] T = U x [vi v 2 v 3 z/ 4 ] T , (3) 

where the mixing matrix is similarly defined to be, 
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U 34 x U 2 4 x U 23 x U u x U 13 x U 12 , 
x U 6 x U t x U t x U t x . 



(4) 



Here the (13) and (14) angles are constrained by chooz to be small: #13, #14 ~ e < eo 
|20|j . The atmospheric neutrino problem constrains #34 (the equivalent of the angle ip m the 
3-flavour case) to be maximal, when both #23, #24 become small. We shall assume they are 
also limited by the same small parameter, eo- 

The other relevant definition is the chosen mass hierarchy in the problem. We define 
the mass squared differences as 5^ = /if — fi 2 , where i,j run over the number of flavours. 
Without loss of generality, we can take 621, £31 (and £41) to be greater than zero; this defines 
the standard hierarchy of masses consistent with the range of the mixing angles, as specified 
above. 



A. Matter effects for neutrinos 



The most important consequence of the highly dense core matter is to cause Mikheyev, 
Smirnov and Wolfenstein (MSW) resonances |2S[] in the neutrino sector. (The chosen mass 



hierarchy prevents such resonances from occurring in the antineutrino sector). In both the 3- 
and 4-flavour cases, this causes the electron neutrino to occur essentially as a pure mass state 
20| , in fact, as the highest mass eigenstate, as can be seen from the schematic illustration 



in Fig. |[ Nonadiabatic transitions near the resonances will alter this result, especially in 
the case of 4-flavours, where there are several level crossings because of the presence of the 
sterile neutrino. Hence we will discuss the purely adiabatic and the non-adiabatic cases 
separately. 



1. The adiabatic case 

We begin with the 3-flavour case. The average transition probability from flavour (3 to 
flavour a is denoted by P a p where a, (3 run over all flavours. It turns out that, independent 
of other mixing angles and the mass squared differences, the survival probability for the 
adiabatic case is, 

Pee = e 2 , (5) 

and hence is small. This is sufficient to determine all the relevant fluxes in the 3-flavour 
case. In the 4-flavour adiabatic case, we again have, 

Pee = e 2 , (6) 

so that the electron neutrino survival probability is flavour independent in the adiabatic 
case. In addition, we also need to know some other probabilities in order to compute the 
fluxes at the detector. We have 

P=P=P=P=e 2 (7) 

1 es 1 se 1 ss 1 ee c V 1 / 
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2. The non-adiabatic case 



Because of the parametrisation it is easy to see that non-adiabatic effects in the form 
of Landau Zener (LZ) jumps are introduced as a result of the values chosen for e and uj. 
The value of e determines whether non-adiabatic jumps are induced at the upper reso- 
nance^) while the value of uj determines whether the non-adiabatic jump occurs at the 
lower resonance. This statement holds both for three and four flavours since in both cases 
the non-adiabaticity in the upper resonance (s) is controlled by e, apart from mass squared 
differences. 

For a large range of e, allowed by the chooz constraint, the evolution of the electron 
neutrino is adiabatic. As a result the lower resonance does not come into the picture at all 
except when e is very close to zero, where the LZ jump probability at the upper resonance(s), 



Ph, abruptly changes to one pT| , |30| . This occurs in a very narrow window: for instance, 
when e changes from 0.02 to 0.01 Ph changes from 0.01 to 0.34 for a mass scale in the range 
of 10 -3 eV 2 . The subsequent discussion is therefore relevant only when e is vanishingly 
small, e <C 10 -2 . (This is not ruled out by the known constraints except lsnd). 

The LZ transition at the lower resonance is determined by the probability, Pl, which is 
a function of the mixing angle uj and the mass squared difference 621- It turns out that Pl 
is zero unless uj is small, sintu ^ 0.2 for mass differences in the solar neutrino range ~ 10~ 5 
eV 2 . This case (of small e, small uj) corresponds to the extreme non-adiabatic limit. When uj 
is large, as in the case of the large-angle MSW or the vacuum solution to the solar neutrino 
problem, non-adiabaticity occurs only at the upper resonance and is therefore partial. 

In our calculations we have used the form for Ph and Pl discussed^] in the appendix of 



Ref. |19| (see also f3"0|). Because of the sharpness of the transition at the upper resonance 
we will set Ph = 1 for small enough e and only consider the dependence of the survival 
probability on the LZ transition at the lower resonance. 

In the 3-flavour non-adiabatic case the only relevant probability is, 

P ee = (l-e 2 )[(l-P L )sin 2 cu + P L cos 2 ^]. (8) 



In the non-adiabatic 4-flavour case we have 



P ee = (l-2e 2 )[(l-P L )sin 2 cc> + P L cos 2 ^] , 

P es = (1 - 2e 2 )[(l - P L ) cos 2 uj + P L sin 2 uj] . (9) 

Since e is small, the flux at the detector is entirely controlled by uj and Pl which is also a 
function of uj. Note that the sum P ee + P es = 1 — 2e 2 is independent of Pl in 4-flavours. Since 
P ee + P es = 1 — P efJ- — P er , this also indicates that the probability of transition of v^ T into v e 
is small in the non-adiabatic case, in contrast to the adiabatic case where P ee + P es = 2e 2 . 



1 The Landau Zener transition probability is defined in terms of the adiabaticity parameter, 7, as 
Plz = ex p[ — ( 7r /2)7P], where F ~ 1. The final expression for 7 given in Appendix B of Ref. [19] 
should be multiplied by the additional factor £31 /60. However, the numerical calculation in that 
paper is correct. 
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B. Matter effects for antineutrinos 

Due to our choice of mass hierarchy, the propagation of antineutrinos is always adiabatic 
(the matter dependent term has the opposite sign as compared to neutrinos). Here also the 
high density of matter in the core causes the electron antineutrino to be produced as a pure 
mass eigenstate, V\. 

In the 3-flavour case, we have 

Pee = (1 - e 2 ) COS 2 LU . (10) 

In the 4-flavour case, we have, 

Pee = (1 -2e 2 )cOS 2 CU , 

P^ = (1 - 2e 2 ) sin 2 to . (11) 

Hence, when u is small, there is very little loss of the V e flux into the sterile channel. Also, 
the sum P^ + Pes = (1 _ 2e 2 ) is similar to the non-adiabatic neutrino propagation so that 
very little v^ T is converted to V e - 

For the NC events, we will also need the following: 

Pse = (1 - 2e 2 ) sin 2 uo + 2e 2 sin 2oo , 

P-= (l-2e 2 )cos 2 u;-2e 2 sm2u; . (12) 
These probabilities can then be used to determine the observed antineutrino fluxes. 



C. The neutrino (antineutrino) fluxes at the detector 

Following Kuo and Pantaleone [[LI]], we denote the flux distribution, d0°/dP, of a neu- 



trino (or antineutrino) of flavour a with energy E produced in the core of the supernova by 
F®. In particular we use the generic label P° for flavours other than v e and V e since 

F2 = F* = ]$ = F° = F? = F«. (13) 

Typically, models for supernovae predict that the v e and v x have hotter spectra than 
v e . This is because the v x decouples first since it has only NC interactions with matter 
and therefore leaves the cooling supernova with the hottest thermal spectrum. The CC 
interactions of u e with matter are larger than those of T7 e and hence u e has the coldest 
spectrum. The model on which this work is based [3Tj] predicts that the average energies 
of the Ue, v e and v x spectra are around 11, 16, and 25 MeV. We make our observations, 
keeping this in mind. 

The v e flux on earth is given in terms of the flux of neutrinos produced in the core of the 
supernova by, 

F = P F° + P F° + P F° 

= PeeFl + (1 " Pee)^ (3-flaVOUrs) , 

= PeeFl + (1 " P ee " P e s)F^ (4-fWirs) , (14) 
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where we have made use of the constraint J2p Pa/3 = 1- This flux is further reduced by an 
overall geometric factor of l/(47r<i 2 ) for the case of a supernova at a distance d from the 
earth. 

Since the probabilities P ee and P es are known, the u e flux can be computed in terms of 
the original fluxes emitted by the supernova. 

The u e flux is independent of the (12) mixing angles u in the adiabatic case. Also, it is 
not very different for 3- and 4-flavours. In both, the observed flux is almost entirely due to 
the original u x flux since P ee = e 2 is small, and is therefore hotter. 

From Eq. ([]) we see that, in contrast, the contribution from the hotter spectrum into 
electron neutrinos in the 4-flavour non-adiabatic case is controlled entirely by e and is small. 
However, the observed flux can be depleted, depending on the value of Pl- The signal in 
the 3-flavour case is drastically different because the contribution of the hotter spectrum 
now depends on e, u and Pl- In general, the possibility of LZ transitions makes the analysis 
more complicated. We will discuss this case numerically later. 

The result for V e is the same, with P ee replaced by Pee, etc. For example, 

F- e = P- e F° + (1 - P^)F° X (3-flavours) , 

= PeeF° + (1 - P^e ~ ^)P° (4-flaVOUrs) . (15) 

There is hardly any mixing of the hotter spectrum into electron antineutrinos in the 4-flavour 
case since (1 — Pee — Pes) = 2e 2 is small. The extent of mixing in 3-flavours depends on the 
value of to through Pee. For small u, there is very little change in the electron antineutrino 
flux (see Eq. ([11])). 

The results are summarised in Tables | and |J for the adiabatic neutrino and antineutrino 
cases, and in Table [ITI] for non-adiabatic neutrino propagation. Expressions for the other 
flavours, needed to compute the NC events, are also given in these tables. For instance, 

2 -F-e ~\~ F T , 

= (P ee + P es + P se + P„)F% + (1 - Pee - P S e)F° e (4 flavours) . (16) 

A similar expression holds for F x with P Q/ g replaced by P-^. 

While the neutral current (NC) combination, (y e + 2v x + V e + 2V X ), remains unaltered 
in 3-flavours, there may be actual loss of spectrum into the sterile channel in the 4-flavour 
case. This should also be a good indicator of the number of flavours involved in the mixing. 

Water as well as heavy water detectors will be sensitive to all these aspects of mixing. 
In the next section, we will discuss the inputs and constraints, both from supernova models 
as well as current neutrino experiments. These will then be used subsequently to predict 
numerically supernova event rates. 



III. INPUTS AND CONSTRAINTS 

We will concentrate mostly on the dominant CC and NC interactions of u e and V e on 
deuteron in heavy water. We will also compare the CC interactions to those at a water 
detector, that is, to T7 e on protons. Interactions on electrons and oxygen nuclei in water and 
heavy water detectors are the same, and have been discussed in detail in Refs. |ITS|,pm. 
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A. Interaction processes and relevant formulae 



Both in water and heavy water detectors, the interactions we are mainly interested in 
are of three types: 

1. Dominant events: These are due to the CC V e p — > e + n interaction in water. In 
heavy water, there are contributions from both V e d — > e + nn and u e d — > e~ pp CC 
processes. At the neutrino energies relevant to this discussion, the cross-sections for 
these processes, and hence the number of events, are about two orders of magnitude 
larger than any of the others. 

2. Electron events: These are due to elastic scatterings of v e , V e , u x , V x on electron targets 
in both water and heavy water. 

3. Oxygen events: These arise from CC scattering of v e and V e on oxygen nuclei in both 
detectors. 

In all these cases, the events are identified by detecting an electron (or positron) in the final 
state. The electron events, although small in number, are highly forward peaked ||32|| . The 
oxygen events have a high energy threshold (E Ue > 15.4 MeV; Ey e > 11.4 MeV) and can 
only occur if there is substantial mixing of the hotter v x spectrum with the v e or the V e 
in which case these events are highly backward peaked p2 |. Both these may therefore be 



readily separated from the mostly isotropic dominant events in water [53]. In the case of 
heavy water, the angular distribution of both the dominant processes is well-known [p T 4]-|36[ 
and is approximately given by 

P{9) oc 1 - -cos# , 
3 

at the energies of interest. Here 9 is the electron (positron) laboratory scattering angle. 
So there are typically twice as many events in the backward direction as in the forward 
direction. This makes the exclusive identification of oxygen events more difficult in a heavy 
water detector. 

A note on the nomenclature "forward" and "backward". We envisage the angular de- 
pendence as being measured in typically six bins of 30° each, the first bin corresponding 
to "forward" and the last one corresponding to the "backward" events. This obviates the 
need for very accurate angle measurements as well as takes into consideration effects like 
electron-rescattering that may smear out the scattering angle. This bin size is also typically 
what is available at SuperKamiokande. 

Since the supernova signal is a short and well-defined signal (lasting about 10 sees), it will 
be possible to detect all these events over background (due to solar and other radioactive 



processes). For instance, in the relevant energy region, SuperKamiokande reports [37] a 
background of 0.2 events/day /kton in a bin of 1 MeV of scattered electron energy from the 
direction of the sun, and half that rate from all other scattering angles. The corresponding 
number for SNO is not available, although the radioactive background is expected to be 



small compared to the solar signal Furthermore, since their angular distribution is 

well-known, angular information on the final state electron (e~ or e + ) may allow us to 
separate out these three types of events. The detailed analysis of signals due to mixing in 



the forward and backward event samples is given in Ref. |20[ for a water Cerenkov detector. 
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These results also hold for a heavy water detector (apart from a scaling factor of 0.9 due to 
the mass difference between water and heavy water). We will not discuss this further and 
concentrate on the new results for the CC events on deuterons. In addition, there are NC 
processes on the deuteron in heavy water that may be detected by the subsequent neutron 
capture and the associated gamma rays. 

The specific processes we will consider, for water and heavy water detectors, therefore 

are 

V e + p — > e + + n , (17) 

and 

T7 e + d—>-e + + n + n , 

v e + d^e~+p + p. (18) 

Common to both detectors are the processes, 

Vol + e — > v a + e [a = e, x, e, x) , (19) 
u e + fO -> e" + X , 

F e + s 6 C -> e + + X , (20) 

where the elastic scattering on electrons involves both CC and NC interactions. 

We will also discuss the interesting possibility of observing NC events on deuteron in a 
heavy water detector: 

V + d — >V + n + p , 

v + d -> v + n + p . (21) 



The cross-sections for Eqs. ([T7|)-(^l]) are well-known p2| , |33| , p2| . The V e p cross-section is large 
in water Cerenkov detectors, being proportional to the square of the antineutrino energy. In 
terms of total number of events, therefore, water Cerenkov detectors are mostly dominated 
by V e p events. The deuteron CC cross sections are comparable though the u e CC reaction on 
deuteron has a slightly lower threshold and a somewhat larger cross-section than the V e CC. 
Also, the V e CC cross-section in heavy water is about 4 times smaller than the corresponding 
one in water due to Pauli suppression. 

When the recoil electron (positron) is detected, the time integrated event rate due to 
neutrinos (or antineutrinos) of flavour a and energy E on target T, as a function of the 
recoil electron (or positron) energy, E e , is as usual given by, 

The index b refers to the time interval within which the (original) thermal neutrino spectrum 
can be assumed to be at a constant temperature Tf,(a). Here tit refers to the number of 
scattering targets (of d, p, e, or O) that are available in the detector. Also, for processes 
involving d and p, the hadron recoil is so small that we assume E e = E — dip, where 5 if 
is the mass difference between the initial and final hadrons (including the binding energy). 



10 



While this results in a small threshold of a few MeV in these cases, the threshold for oxygen 
processes is greater than 10 MeV. The total number of events from a given flavour of neutrino 
in a given bin, k, of electron energy (which we choose to be of width 1 MeV) then is 

rk+l A]V T 

w-L AE -Ji- (23) 

In the case of NC events in heavy water, the total number of events is calculated according 

to 



£ At b JdEF a (b)a NC (E) . (24) 



Here again the cross-section is well-known [33]. 



B. The supernova flux inputs 



As in [Jl9| , p0"| , we compute the time integrated event rate at prototype 1 kton water 
and heavy water detectors from neutrinos emitted by a supernova exploding 10 kpc away. 
Results for any other supernova explosion may be obtained by scaling the event rate by the 
appropriate distance to the supernova and the size of the detector, as shown in Eq. (p2|). 
We assume the efficiency and resolution of the detectors to be perfect; this will only slightly 



enhance the event rates near the detector threshold [20]. 

We use the luminosity and average energy distributions (as functions of time) for neutri- 
nos of flavour a and energy E as given in Totani et al. ||31||, based on the numerical modelling 



of Mayle, Wilson and Schramm ||39|| . In a short time interval, At b , the temperature can be 
set to a constant, T b (a). Then the neutrino number flux can be described, in this time 
interval, by a thermal Fermi Dirac distribution, 

f °'">^°|^ (e Xp ( g/ £ a )) + 1 ) - < 25 > 

at a time t after the core bounce. Here b refers to the time-bin, t = t + bAt. We set 
the time of bounce, to — 0. The overall normalisation, N Q , is fixed by requiring that the 
total energy emitted per unit time equals the luminosity, £ b (a), in that time interval. The 
precise values for £ b (a) and T b (a) are taken from Ref. pH] . The total emitted energy in 
all flavours of neutrinos is about 2.7 x 10 53 ergs. The general features of the model are 
as follows: The temperature is roughly constant over the entire period of emission (lasting 
roughly 10 seconds). Typical values are T b (a) = 3.15(E) a , with the average energy of each 
flavour, (E) a = 11,16,25 MeV for a = v e ,V e and u x respectively. Also, the total emitted 
energy is more or less equally distributed in all flavours. (The luminosity of v e is higher than 
that of other flavours at early times, while that of v x is higher after 1 sec). The number 
of neutrinos emitted in each flavour, however, is not the same since their average energies 
are different. While we use the values for the time dependent temperature and luminosity 
as given by Ref. PT| in our analysis, we also examine the effects due to possible variations 



of these parameters and hence the sensitivity of our results to the details of the supernova 
model used. 
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With large matter effects present in both the neutrino and antineutrino sector, the va- 
lidity of the average energies of v e , V e and u x as 11, 16 and 25 MeV, respectively, which were 
calculated without mixing, may be questioned. This is especially so because re-scattering 
effects involve the flavour states which may then equilibrate at different temperatures. We 
however note that the highly dense matter projects the u e and V e states as almost pure 
mass eigenstates. Hence thermalisation is not affected by effects of mixing. In fact, the ef- 
fects of mixing are significant only when the resonant densities are reached, when the MSW 
effect can mix different flavour states. For the parameter values as allowed from current 
neutrino data, this occurs only outside the neutrinosphere (R ~ 10 4 km), and not at the 
core (R ~ 50 km) where most of the neutrinos are produced. This mixing therefore occurs 
between spectra which are already thermalised with the above-mentioned temperatures. In 
the case of and u T this argument does not go through. In the highly dense core, these are 
mixtures of more than one mass eigenstate. However, both mix only into each other and 
scatter through exactly the same processes. Hence their temperatures also remain the same 
as in the no-mixing case. 



C. The mixing parameters 

We impose the following known constraints on the mixing matrix in vacuum both for 
three and four flavour scenarios. Consistent with the chooz constraint, namely sin0 ~ 
e < e = 0.16 which we have imposed at the level of the parametrisation itself, we choose 
e = 0.08 for the 3-flavour adiabatic case and e = 10~ 4 ~ for the non-adiabatic case. In the 
case of the 4-flavour scheme, we set 6*13, 6*14, #23, #24 = £■ 

The constraint from the atmospheric neutrino analysis implies that the relevant angle 
ip{= #34) ~ vr/4, is near maximal and the relevant mass squared difference is of the order of 
10 -3 eV 2 . Neither of these constraints directly enter our calculations except to determine 
whether the upper resonance is adiabatic or not depending on the value of e as constrained 
by the chooz findings. We consider both possibilities here. 

The value of the (12) mixing angle, lu, is not yet known. Combined data on solar 
neutrinos give three possible values [|4(J : 

1. sm 2 2u = 6.0 x 10- 3 ,5 2 i = 5.4 x 10" 6 eV 2 (SMA). The small angle MSW solution. 

2. sm 2 2u = 0.76, 5 2 i = 1.8 x 10~ 5 eV 2 (LMA). The large angle MSW solution. 

3. sin 2 2cj = 0.96,^21 = 7.9 x 10~ 8 eV 2 (LMA-V). The large angle vacuum solution. 



This has been slightly modified p4| , f42|| in view of new data from SuperK; however, it remains 
true, in general, that u may be small or large, with 621 < 10~ 5 eV 2 . 
We will now discuss the results numerically for all these choices. 



IV. NUMERICAL RESULTS 

The numerical calculations are done by following the evolution of the mass eigenstates 
through all the resonances including the appropriate jump probabilities when the transition 
is non-adiabatic (when e is small). Furthermore, the LZ jump at the lower resonances is 
significant only for small values of u. 
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A. The electron (positron) spectrum 



The CC event rates on deuteron computed using the above inputs are displayed in Figs. 
H The time-integrated event rates per unit electron energy bin (of 1 MeV) are shown as 
a function of the energy E e of the detected electron. The solid lines in all the figures refer 
to the case when there is no mixing and serve as a reference. Results for 3- and 4-flavour 
mixing are displayed in each of these figures as dotted and dashed lines respectively. 

Fig. |3] shows the predictions for the (u- independent) adiabatic v e d CC interaction when 
e = 0.08. Mixing enhances the high E e event rates for both 3- and 4-flavour mixing, which 
cannot be distinguished here. Furthermore, mixing shifts the peak of the spectrum to higher 
energies, from 15 MeV to about 28 MeV because of the admixture of the hotter v x spectrum 
and its significantly higher cross-section with deuterium. This high energy shift should be 
clearly observable. 

Nonadiabaticity at the upper resonance occurs when e ~ 0. Then the adiabaticity at 
the lower resonance is determined by the value of uj. In fact, Pl — unless uj is small, 
sin a; ~ 0.2. Hence this is relevant only for the small angle MSW solution (SMA). We show 
the dependence of Pl on the neutrino energy for the small-cu 3-flavour case in Fig. [|. The 
4-flavour result is similar, with a scale factor roughly 0.7. It is seen that Pl increases with 
energyalthough it still does not reach unity for the relevant supernova neutrino energies. 
For larger uj, that is, for the large angle MSW (LMA) and the vacuum (LMA-V) solutions, 
Pl = 0. 

In Fig. |5] we show the electron spectrum for the non-adiabatic v e d CC interaction when 
e is small, in fact near zero. The fully non-adiabatic case corresponding to small uj is shown 
in comparison with that for a larger value of uj in the figure. Here, 3- and 4-flavour mixing 
give drastically different results, but the small-e scenario is in general not very sensitive to 
the chosen values of uj. 

In Fig. ^ we show the positron spectrum due to V e d CC interactions for two different 
choices of uj. (Here e = 0.08, but the results are essentially the same even if e is nearly 
zero since this sector is always adiabatic.) Mixing has appreciable effects only for large uj; 
however, mixing does not affect the peak position, unlike in the adiabatic v e case. 

We would like to point out that the pure v e spectrum may not be observable in heavy 
water. It may be separated out from the total events sample if the V e spectrum can be 
reliably separated out by various detection techniques such as looking for two neutrons in 
coincidence with the positronf]. However, we will show that the total number of events will 
still be sensitive to the mixing parameters. We will therefore discuss both the total as well 
as the individual v e and v e spectra. 

The total event rates due to the dominant CC v e d and V e d processes are shown in Figs. [7| 
and H for the two choices of e. The total isotropic CC events in water, due to T7 e p alone, 
are also shown for comparison. The most significant difference between the two is that of 
the adiabatic case with small uj (the lower two panels of Fig. 0) which is independent of the 
number of flavours. This is because of the enhancement in the v e events, independent of uj. 
At all u, the peak is at a higher energy than expected from the no-mixing case in heavy 



2 We thank the referee for pointing out this possibility. 



13 



water but remains the same for a water detector. This shift may be sufficiently significant 
and therefore observable, in the adiabatic scenario, for all u. Finally, the upper two panels 
of Fig |7| indicate that a significant depletion in the observed events in water, together with 
an enhanced number of events in heavy water is an unambiguous signal of 4-flavour mixing 
with large u (dashed lines in Fig. [?[). 

The corresponding results for the dominant CC events in the non-adiabatic case are 
shown in Fig. ^. Here there is no significant shift in the spectral peak. Also, the signals in 
water and heavy water are very similar, with the signals being either enhanced or depleted 
similarly in both. For instance, the large us 4-flavour signal shows depletion both in water 
and heavy water. This may be difficult to distinguish from the no-mixing case if the overall 
normalisation of the supernova spectrum is uncertain by more than a factor of two. In all 
cases, the small-e, small-tu scenario also cannot be distinguished from the no-mixing case. 
Hence the small-e case may be difficult to establish unambiguously, independently of the 
supernova model inputs. 

Keeping in mind that the supernova dynamics may have large uncertainties, we will later 
also analyse the ratios of the total number of events in water and heavy water detectors. 
These are likely to be less sensitive to the supernova models (although they do depend on the 
temperature hierarchy for v e , V e and u x ) and hence may be more robust signals of mixing. 



B. Integrated number of events 

The predicted time integrated number of events resulting in a scattered electron with 
energy, E e > 5 MeV (which is a typical threshold for Cerenkov detectors), are shown in 
Tables [IV| and [V|, for the adiabatic and non-adiabatic cases respectively. As before, the 



number is calculated assuming a supernova explosion at 10 kpc for a 1 kton detector. Listed 
are the dominant CC events on deuterons in heavy water and free protons in water, along 
with the elastic scattering events on electrons and the CC events on oxygen nuclei. 

For heavy water, we have listed the individual contributions from u e and V e on deuteron. 
In water, the corresponding dominant events are from V e on p. The elastic events are from 
v e e, V e e, v x e, and V x e. Since they will all be detected in the extreme forward direction, they 
have been summed up and listed as total v e events in the tables. The oxygen events, listed 
as v O, include both v e O and F e O CC events, which will predominantly be in the backward 
direction, especially when enhanced by mixing. In particular we tabulate the events for 3- 
and 4-flavour mixing when uj is both small and large. 



In Table [TV], we show the results for the fully adiabatic case when e = 0.08. It is 
seen that the bulk of the events (more than 90%) are from the dominant events on p or 
d. Mixing always enhances the v e d channel by more than a factor of two; hence adiabatic 
propagation always predicts an enhanced rate of total events in heavy water even though 
there is reduction in the V e d channel (the other dominant process) for some parameters. 
In contrast, the total number of events in water may even go down as compared to the 
no-mixing case, depending on the parameter values. 

In Table [V], we show similar predictions for the case when e ~ or when the upper 
resonance(s) become fully non-adiabatic. The value of u then determines whether or not the 
lower resonance is adiabatic. Recall that the antineutrino propagation is always adiabatic. 
Again, the contribution from the ue and vO events is small compared to the dominant 
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events. Small changes in the V e p and V e d events between Tables |V| and [V| are due to 
changes in the value of e. Irrespective of the parameter values, it is seen that there is never 
any depletion in the 3-flavour case. As in Fig. [8] an interesting scenario occurs when to is 
large in the 4-flavour case. Here the total number of dominant events in both v e d and V e d 
in heavy water and in V e p in water, is reduced by a factor proportional to cos 2 a;. This is 
the only scenario where there is depletion in both water and heavy water. 



C. Possible discrimination of various "mixing models" 

So far, we have assumed that most of the mixing parameters are known and used the su- 
pernova measurement as a potential check for self-consistency of the model parameters. This 
is because there is still very little known about the supernova neutrino spectrum through 
observations and hence there is both theoretical and experimental uncertainty about the 
details of the neutrino spectrum. However, it is still instructive to actually turn the ques- 
tion around and ask, suppose another supernova explosion is observed through its neutrino 
emission. Will an excess over (or a depletion from) the expected number of events unam- 
biguously determine some of the model mixing parameters ? The answer to this question 



can be obtained from Table |VT| . It is seen that certain classes of models may be ruled out, 
depending on the observation. We define the ratio of the total number of events (from all 
possible interactions with an electron (or positron) in the final state) potentially observed 
from a future supernova (equivalently, the prediction from a given neutrino mixing model) 
to the expected number of events without mixing (for E e > 5 MeV): 

observed number of events . . 

1 calculated number without mixing ' 

where i = D, H refer to 1 kton heavy water and water detectors respectively. The denomi- 
nator refers to the the expected number of events (using a standard supernova model, with 
no mixing) as computed from a Monte Carlo simulation that takes into account detector 
resolution, efficiency, etc., consistent with the detector at which the events were observed. 
An observation may find Ri > 1, Ri < 1 or Ri ~ 1. Note that even though the ratio refers to 
the total number of events, the inferences drawn reflect mainly the behaviour expected from 
the dominant CC processes on protons in water and deuterons in heavy water. The mixing 
models (with model parameters e and u, including adiabaticity) consistent with, or predict- 
ing, such an observation are shown in Table [VT]. Here the non-adiabaticity, that is, the value 
of Pl at the lower resonance has been computed assuming a typical value of 5 2 i = 10~ 5 eV 2 . 
Water detectors cannot distinguish adiabatic (A) and non-adiabatic (N) scenarios, that is, 
whether or not e is different from zero, but can distinguish the number of flavours when u is 
large (see the last column of Table |VTD . In D 2 0, however, most models predict Ri > 1. The 



only observation of Ri < 1 in D 2 occurs for the 4-flavour non-adiabatic case with e ~ 
and large u. 

On combining data from water and heavy water detectors, an improved discrimination 
of model parameters is possible, as can be seen from Table |VT|. Here the different values 
of Rh and Rjj are listed, along with the models that are consistent with such a combined 
observation. First of all, it is seen that combining the two measurements immediately 
allows for a separation between the adiabatic and non-adiabatic cases and hence whether 
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e is different from zero, except in the 3-flavour case with large uj. It should be noted that 
the e, uj — > scenario is unlikely to be distinguished from the no-mixing case. It is also seen 
that certain combinations of Rd and Rh do not occur for any of the allowed parameters 
values. For instance a depletion in heavy water is only consistent with depletion in water 
as well. Any other result in water indicates that the overall normalisation of the spectrum 
is probably in error. An occurrence of such "forbidden" combinations may therefore can 
be used as a check on the overall normalisation of the supernova spectrum. This result of 
course is limited to the class of models we are analysing here. 

The following scenario is best suited to determining the value of uj. (1) There are fewer 
isotropic events than expected (Rh < 1) in a water Cerenkov detector such as SuperK. This 
reduction factor determines cos 2 u. (2) The same reduction factor (Rd < 1) fits the data 
from a heavy water detector such as SNO. This can imply that the correct mixing matrix 
is one with 4 flavours and large uj, with non-adiabatic neutrino propagation. (3) If on the 
other hand there are enhanced number of events (Rd > 1) at the heavy water detector, it 
clearly indicates adiabatic neutrino propagation. This in turn implies that e is different from 
zero which has so far been claimed only by lsnd. Variation of supernova input parameters, 
which we will discuss in the next section, does not alter this result. 

These qualitative features can be quantified by defining the double ratio, 

Rd/h = , (27) 

which is independent of the overall normalisation of the neutrino flux and hence provides a 
better diagnostic. In practice, it may not be possible to directly take a ratio of the data from 
water and heavy water detectors since the two measurements will differ in their systematics, 
apart from such considerations as detector efficiency and resolution. Since Ri,i = D,H are 
normalised to the theoretical expectancy including these considerations, Rd/h is not likely 
to be sensitive to details of detector design and can thus provide a robust, quantitative 
indicator of different types of mixing. 

This double ratio (where Ri has been calculated as before) has been shown in Figs. |9| 



and [10] as a function of uj. The ratio is plotted for the total number of events with the cut 
on the observed electron (positron) energy, E e > 5 MeV in Fig. |9|. The two curves in each 
figure correspond to the two values of e when the propagation at the upper resonance (s) 
is purely adiabatic, e = 0.08 (solid lines) or purely non-adiabatic, when e = 10~ 4 (dashed 
lines). These are the two cases that a water detector is normally not able to resolve. Also 
shown are dotted vertical lines corresponding to the solutions allowed by the solar neutrino 
problem: sin 2 2uj = 0.006,0.76 and 0.96. Non-adiabaticity at the lower resonance has been 
computed using 8 2 \ = 10~ 5 eV 2 as before. While the 3-flavour mixing case is shown on the 
left, the 4-flavour result is plotted on the right. 

Obviously, a value of unity is expected for the case of no-mixing. We analyse each case 
in turn. 

1. We see from Fig. |9| that the double ratio Rd/h is always strictly greater than one for 
the adiabatic case, independent of uj or the number of flavours, /. 

2. Even in the non-adiabatic case, it can be less than one only when / = 4. Note however 
that currently allowed values of uj lead to Rd/h ~ 1 in the non-adiabatic case. The 
case Rd/h < 1 occurs only for intermediate values of uj. 
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3. For / = 3 (left panel of Fig. ^j, the double ratio at small u is different for the adiabatic 
and non-adiabatic cases, which may therefore be distinguished. However, it may be 
very difficult to distinguish these two for large values of u. 

4. On the other hand, these two cases are easily distinguished for all u for / = 4 as can 
be seen from the right panel in the figure. 

5. While the double ratio is similar for both / = 3, 4 for small u, sin 2 2uj < 0.1, the number 
of flavours can be distinguished for larger values of u, especially in the currently allowed 
region, only in the adiabatic case. 

6. However, as stated before, independent of /, the small u non-adiabatic solution with 
e ~ cannot be distinguished from the no-mixing case. 

Keeping in mind that a thermal neutrino flux distribution such as the one we have used may 
overestimate the high energy spectrum, the case 5 < E e (MeV) < 40 is shown separately 
in Fig. [l(| Most of the features survive the cuts; hence this ratio is likely to be a stable 
indicator of mixing. 

Finally, we note that the denominator of the double ratio is dominated by V e events. 
Hence, the same discriminatory power can be achieved using data from a heavy water 
detector alone if it is possible to separate the v e d and V e d events. As stated earlier, this 
may be possible, for instance, at SNO, by detecting both the neutrons in coincidence with 
the positron emitted in the V e d interaction. SNO is also planning to increase the neutron 



detection efficiency (to more than 80%) by adding salt to the heavy water fl38[. In this case, 
the double ratio, defined for heavy water alone, 

U = %^ , (28) 



will provide as much information as the double ratio Rd/h- Here, 

Re 



observed number of v e events 



calculated number without mixing 

while Re+e is a similar ratio, defined for the total number of events from both v e and V e 
interaction with deuterium: 



R P 



observed number of v e and v e events 
calculated number without mixing 



Both R e and R e +e are calculated for a heavy water detector. To a very good approximation, 
we have 

Rd/h « r d , (29) 

The approximation arises partly from ignoring events due to electron and oxygen targets 
in rd- The error also arises from the differences in the denominators of the two ratios, one 
involving V e p and the other T7 e d. Despite a mild energy dependence of the ratio of these 



two cross-sections [33,35], it turns out that the ratio of the total events expected from these 



two processes remains in the range 4 ± 0.1 (see Table |TVD . This is true both when there 
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is no mixing, and with mixing, for all allowed values of mixing angles (<f),w). Hence this 
factor cancels when the ratio is expressed in terms of Rh- Hence the approximation in 
Eq. p9| should be valid to within a few percent. In addition, the double ratio will have the 
advantage of reduced systematic errors, since data from different experiments do not have 
to be combined in order to calculate it. Hence it will be useful to calculate such a ratio, by 
separating out the V e events on deuteron. 



D. Sensitivity to supernova model parameters 

So far, we have discussed the sensitivity of the supernova neutrino spectrum to various 
neutrino mixing parameters. However, the supernova model parameters (temperature and 
luminosity) are themselves uncertain and still need to be experimentally established. It is 
therefore important to study the effect of variation of these parameters on the results we 
have so far obtained. 

Supernova dynamics is a very complicated issue. Here we will follow a simple-minded 
approach. Changes in the luminosity affect the overall normalisation while changes in the 
temperature (or average energy) change in the shape of the spectrum. Variations in these 
parameters, while being time-dependent, are not random, but systematic. For instance, the 
supernova model parameters depend on the protoneutron star mass (an increase of which 
increases both the average energy and luminosity of neutrinos) as well as the underlying 
high-density equation of state and the initial conditions. The effect of this on the total 
neutrino spectrum has been studied in Ref. f43fl . (The temperature variation in the spectra 
of individual flavours is not discussed). The study indicates that the typical temperature 
variation of the total neutrino spectrum at all times does not exceed about ±1 MeV. While 
the variation due to uncertainties in the initial conditions is relatively small, the average 
energy systematically decreases for smaller protoneutron mass stars and those evolving with 
a stiffer equation of state. On the other hand, the luminosities are virtually identical for all 
these cases until a time t ~ 10 sees., by which time most of the detectable neutrinos are 
emitted. 

We will first estimate the systematic errors due to uncertainties in the supernova tem- 
perature. In the absence of detailed information on temperature variation of the individual 
flavours, we shall assume a systematic time-independent increase (or decrease) in the tem- 
perature of spectra of all flavours by 1 MeV. This will then be an estimator of the outer 
limits of variation of the results from the original calculation. 



Fig. |TT| shows the expected number of events due to v e d interaction in the absence of 
mixing and when the temperature is systematically increased or decreased by 1 MeV in all 
time bins, 

T b {a) -> T b (a) ± 1 MeV . 

The base-line supernova spectrum is shown in comparison as a solid line. There is of course 
a shift in the spectral peak (by around 3 or 4 MeV); however, there is a large change in 
the high energy part of the spectrum (accentuated due to the E 2 dependence of the cross- 
section). It will still be possible to distinguish the adiabatic mixing case since the increase at 
high energy in this case is substantially larger than from errors in the supernova spectrum. 
However, other cases, especially the non-adiabatic cases, will not be clearly distinguishable. 
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It must be noted that in any event the spectral peak for v e events is a good index of the 
temperature of the spectrum, either of the unmixed v e or of the hot v x spectrum. 



Fig. [12] shows the results for the case of the V e d unmixed spectrum. There is a similar 
dependence (since the energy dependence of the cross-section is the same as in the v e d case. 
Since mixing does not significantly shift the V e spectral peak (see Fig. ||), this will remain a 
good indicator of the corresponding spectral temperature. 

While there are large variations in the results for the individual spectra, these will be 
cancelled out in the double ratio Rd/h of the events in heavy water and water (in fact, 
this is the purpose of constructing such a ratio). This can be seen from Fig. [T3], where the 
same double ratio, Rd/h, is plotted for the different temperature sets, T (as solid lines) and 
Til (as dashed lines). It is likely that any time-dependence of the temperature variation 
that we have ignored, will affect the numerator and denominator of the ratio in the same 
way; hence inclusion of time dependence should not affect this analysis. In computing this 
ratio, the "observed number of events" as required for the calculation of Ri, i = D,H, 
in Eq. ^6] is now determined both by the mixing parameters as well as by the modified 
supernova model. It is seen that there is very little sensitivity to the variations in the 
temperature. This is especially so in the adiabatic case. A greater sensitivity to the model 
parameters in the small u, small e (non-adiabatic) case occurs because of the presence of 
the additional energy- dependent factor, the Landau-Zener transition probability, P^. Hence 
inclusion of temperature variations in the supernova model does not change the conclusions 
about discrimination of different mixing models. 

We add a note on variations in the luminosity due to, for example, choice of different 



initial conditions [|3|. This affects the overall normalisation which is an indicator of the 



total energy emitted. However, the double ratio will not be sensitive to this, unless these 
changes are extremely time-dependent. It is doubtful whether reasonable conclusions can be 
drawn in such a case, unless there are significantly large numbers of events in each time bin. 
In short, while the individual flavour spectra can be significantly modified by uncertainties 
in the supernova model parameters, the double ratio Rd/h is largely insensitive to such 
variations. Hence it is a good indicator of mixing. 

A final remark about statistical errors. As already stated, the background to a supernova 
signal due to both radioactivity as well as solar events is small at SuperKamiokande and 
SNO. Hence the signal will be clearly defined. In any case, the statistical errors (assuming a 
1 / \/N error for both the numerator and the denominator and adding suitably in quadrature) 
have been calculated for the double ratio Rd/h- This has also been shown in Fig. [T3|. The 
errors are so small that they are not visible, except as a slight thickening of the lines near 
the oj — > 1 region. Of course, if the supernova is 50 kpc and not 10 kpc away, the statistical 
errors become 5 times larger. Recall however that we have computed the events in 1 kton of 
detector. Larger detector volumes will further reduce this error. In general, the statistical 
quality of the signal, while being good, will depend on the size of the detector as well as the 
distance to the supernova. 

E. Neutral Current events 

As is well-known, heavy water detectors can directly observe NC events. This is very 
important in the context of supernova neutrinos since neutrino emission from supernovae 
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is practically the only observable system where neutrinos (and antineutrinos) of all flavours 
are emitted in roughly equal proportions. Note that there are also NC events on oxygen 
targets in both water and heavy water, with a characteristic signal of photons with energies 
in the range of 5-10 MeV ||44|| . However, these events are fewer in number than the NC 
events on the deuteron that we will discuss here. 

While there is no loss of NC events in the case of 3-flavour mixing, the existence of 
a fourth flavour will be signalled by loss of NC events into this sterile channel. This can 
be seen from Tables |VIII| and OCl where the total number of NC events from neutrinos or 



antineutrinos of all flavours (with E v > 3 MeV) is listed for different possible values of u> 
consistent with the solar neutrino expectation. It is seen that the number of NC events is 
not very sensitive to the value of u; however, from Table [Villi it is clear that in the adiabatic 
case, there is about 25% depletion with 4-flavour mixing when compared to the no-mixing 
case. If the value of e and the overall normalisation of the spectrum is known, NC current 
events can be used to discriminate between three and four flavour mixing. Recall that the 
CC events are always enhanced by a factor of 1.5-2 for the adiabatic case. Hence more 
conservatively, the NC events can be used to normalise the supernova spectrum to at least 
within 25%. 

When the upper resonance is non-adiabatic, however, part of the signal gets regenerated, 
especially for smaller u. Hence, as Table [X| shows, there are roughly the same number of 
NC events with and without mixing in the non-adiabatic case, independent of the number 
of flavours. The corresponding CC channel shows severe depletion only when u is large; 
otherwise, it is either enhanced, or the same as the no-mixing case. Hence here also the NC 
events can be used to normalise the supernova spectrum. 



V. 4-FLAVOUR (3 + 1) MIXING SCHEME 

So far, all results in the 4-flavour analysis referred to the (2 + 2) scheme as shown in 
Fig. [T]b. We briefly discuss results in the (3 + 1) scheme where the mixing matrix is defined 
through: 

We ^ v T v s ] T = U x [v x v 2 u 3 z/ 4 ] T , (30) 

and 

U = U u x C/ 24 x U 2 3 x U u x U 13 x U12 , 

= U p xU € xU lp xU e xU e xU u} . (31) 

Here the (13) and (14) angles are constrained by chooz to be small: #13, 9 14 ~ e < eo [0 as 
in the (2 + 2) scheme. The atmospheric neutrino problem now constrains #23 (the equivalent 
of the angle if) in the 3-flavour case) to be maximal, when # 24 becomes small, #24 < eo- 
However, the (34) mixing angle p is not constrained by any known experimental data. 

Since u e is produced in the supernova core in essentially the V4 mass eigenstate, any non- 
adiabaticity results in jumps near the upper MSW resonances. The adiabaticity parameter 
here will be determined by the #14 angle, which is again small, # 14 = e. However, the 
adiabaticity parameter at the lower resonance depends on the unknown angle p and hence 
we do not comment on the non-adiabatic case here. 
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The expressions for the fluxes as observed on earth are given in Eqs. (|I4]) , ([15]) and ([16]) 
with probabilities P a p computed for the (3 + 1) scheme. The relevant probabilities, P ee and 
P es , for the CC events in the adiabatic case are independent of the unknown angle p and are 
in fact the same in both the (2 + 2) and (3 + 1) schemes. Hence all the results (as shown 
in Figs. |3], |6], |7|) for the CC events on deuteron in the adiabatic sector are insensitive to the 
position of the sterile neutrino in the case of 4-flavours. This is true for the CC events on 
oxygen as well. 

In the NC case, we need the probabilities P se and P ss which are different from the (2 + 2) 
case. Ignoring small terms of order 0(e), we have, 

4' 

44/2 , 

44/2 , (32) 

where c p , s p are cos p and sin p respectively. 

While these fluxes do depend on p it turns out that the suppression factor in the adiabatic 
case is again around 75% for both large and small values of u, again independent of the 
value of p. This is because the dominant contribution to the NC sector is from the F x and F% 
fluxes, as in the (2 + 2) case. These terms are very weakly dependent on the unknown angle 
p. Hence in the NC sector as well, the (adiabatic) (3 + 1) 4-flavour scheme gives almost the 
same predictions as the (2 + 2) scheme. Because of this, there will not be much difference 
in the elastic events on electrons as well. 



p 

P— = 

1 se 

P— = 

1 ss 



VI. SUMMARY AND DISCUSSION 

To summarise, we have contrasted signals from supernova neutrinos (antineutrinos) in 
water and heavy water detectors. We include the dominant charged current events from 
deuteron targets in heavy water, and proton targets in water, as well as the elastic scattering 
off electrons and charged current events on oxygen in both detectors. In all cases, an electron 
(or positron) is detected in the final state. The detailed distribution of events as a function 
of the scattered electron (positron) energy depends on the number of flavours and the mixing 
parameters in a complex manner. We have discussed all these cases. 



In particular, Figs. |9] and |10] show the combined sensitivity of water and heavy water 
detectors to the neutrino mixing parameters for such events by defining a double ratio of 
the observed to expected number of events in heavy water and water detectors respectively. 
These results reflect essentially the behaviour of the dominant charged current v e and V e 
events on deuterons which are comparable to the V e charged current interaction on protons 
in water. However, its dependence on the mixing parameters is very different from that for 
water. It turns out that a comparison of the signals from water and heavy water detectors 
can yield important information on not only mixing parameters but also on the number 
of flavours involved. While 3-flavour mixing typically results in an enhanced event rate, 
4-flavour mixing can lead to substantial decrease in the number of events, depending on the 
mixing parameters. 
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Furthermore, we have performed a simple-minded analysis of the systematic errors in- 
volved due to uncertainties in the supernova model parameters. We have shown that the 
double ratio that quantifies the relative variation due to mixing in water and heavy water 
detectors is largely insensitive to variations in the supernova model parameters (temperature 
and luminosity) used. 

We have also briefly discussed the neutral current events in heavy water. These signals 
may facilitate determination of the overall normalisation of the supernova neutrino spectra. 
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TABLES 



No. of flavours 


Neutrino flux at detector, Ff 


3 
4 
3 
4 


F e = e 2 F° + (1 - 

F e = e 2 F e ° + (1 - 2e 2 )F° 

2F X = (1 + e 2 )F° + (1 - e 2 )F e ° 

2F X = (4e 2 )F° + (1 - 2e 2 )F e ° 


TABLE I. v e , v x = v^v T fluxes at the detector in the extreme adiabatic limit for 3- and 
4-flavour mixing. 


No. of flavours 


Antineutrino flux at detector, Ff 


3 
4 
3 
4 


F S = (1- e 2 )c 2 Fl + (4 + e 2 c 2 )F° 
F g = (l-2e 2 )c 2 F° + (2 e 2 )F° 
2F- X = (1 + cl - e 2 c 2 )F° + (4 + e 2 c 2 )F° 
2F 5 = (2 - 4e 2 )F° + +2e 2 (l - s 2u; )F| 


TABLE II. v e , v x = u^, v T fluxes at the detector for 3- and 4-flavour mixing. Here c w = cos a;, 
= sinco. 


No. of flavours 


Neutrino flux at detector, Ff 


3 

4 
3 

4 


F e = (l-e 2 )[(l-P L )sl + P L cl]F2 
+[ l_ ( l_ e 2 )((1 _p L)s 2 + p LC 2 )]F o 

F e = (1 - 26 2 )[(1 - P L )sl + P L c 2 ]F e ° + 2e 2 F° 
2F X = [1 + (1 - e 2 )((l - P L ) S 2 + Pl4)]F° 

+ [ l_(l_ e 2 )((1 _p L)s 2 + p LC 2 )]jP 

2F X = 2(1 - 2e 2 )F° + 2e 2 [1 + (1 - 2P L )s 2u) ] F e ° 



TABLE III. Neutrino fluxes at the detector when non-adiabatic effects are introduced. While 



the transition is assumed to be fully non-adiabatic at the upper resonances, it is controlled by the 
jump probability Pl at the lower resonance. Here c w , and S2u> refer to cos a;, sinw, and sin2u; 
respectively. 
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Heavy water 


Water 




l\IO 


3-flavours, s\ u = 


4-flavours, 


„2 

s 2cj — 


J.NO 


3-navours, 


„2 

b 2u — 


4-flavours, 


„2 

6 2oj — 




Mixing 


0.96 


0.006 


0.96 


0.006 


Mixing 


0.96 


0.006 


0.96 


0.006 


V e d(p) 


72 


183 


183 


181 


181 

















V e d{p) 


71 


85 


71 


43 


71 


290 


329 


291 


177 


291 


v e 


8 


9 


9 


8 


8 


9 


10 


10 


9 


9 


v 


5 


29 


26 


24 


26 


5 


32 


28 


27 


29 


Total 


156 


306 


289 


256 


286 


304 


371 


329 


213 


329 



TABLE IV. Total number of events in 1 kton D2O with electrons (positrons) in the final state 
with energy, E e > 5 MeV. Listed are the contributions from CC events on deuterons due to both v e 
and u e , events due to the elastic scattering of all flavours and antiflavours of neutrinos on electrons 
(labelled ve), and the CC events from v e and V e scattering on oxygen nuclei (labelled uO). The 
results due to no-mixing, and mixing with 3- and 4-flavours in the adiabatic case with e = 0.08 are 
shown in the three columns. The results with 3- and 4-flavour mixing are shown for two values of 
lu: u large (sin 2 2uj = 0.96) and to small (sin 2 2a; = 0.006). For comparison, the total number of 
events in water are also listed for the same set of model mixing parameters. The deuteron target 
is replaced by free protons here. 





Heavy water 


Water 




No 


3 flavours, 


b 2u — 


4 flavours, 


s 2 - 


No 


3 flavours 


s 2 - 


4 flavours, 


b 2u) — 




Mixing 


0.96 


0.006 


0.96 


0.006 


Mixing 


0.96 


0.006 


0.96 


0.006 


v e d{p) 


72 


139 


88 


29 


52 

















V e d(p) 


71 


85 


71 


43 


71 


290 


329 


290 


174 


290 


v e 


8 


9 


8 


5 


7 


9 


10 


9 


6 


8 


v O 


5 


21 


8 


3 


4 


5 


23 


8 


3 


5 


Total 


156 


254 


175 


80 


134 


304 


362 


307 


183 


303 



TABLE V. The same as Table [V but for a small value of e 
nance(s) is non-adiabatic. 



10 4 so that the upper reso- 
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Models allowed by the corres 
D 2 


ponding value of R measured in i 
H 2 


R > 1 


(3uj l ) a ,n, (4wl)a, (S,4ws)a 


(3ll>l)a,n 


R < 1 




(4ul)a,n 


R ~ 1 


No mix, (3, 4:Ll>s)n 


No mix, (3, Alo s )a,n 



TABLE VI. List of neutrino models which can be discriminated by the D 2 and H 2 detectors 
from different values of R. Here R is defined as the ratio of the observed number of events with 
electrons of energy E e > 5 MeV in the final state, to the calculated number without mixing. 
The various models are specified by the number of flavours 3 or 4, by the value of uj (ujl and 
los refer to sin 2 2uj = 0.96,0.006 respectively), and by the suffix A and N referring to adiabatic 
and non-adiabatic propagation at the upper resonance(s), corresponding to e much larger or much 
smaller than 10 -2 respectively. 





Models which are allowed 


Rd > 1; 


> 1 


(3ll>l)a,n 


Rd > 1; 


An < 1 




Rd > 1; 


ite ~ 1 


(3,^s)a 


Rd < 1; 


Rh > 1 


None 


Rd < 1; 


< 1 


(4^l)at 


Rd < 1; 


R H ~ 1 


None 


Rd ~ 1; 


> 1 


None 


Rd ~ 1; 


< 1 


None 


Rd ~ 1; 


~ 1 


No mixing, (3,Aujs)n 



TABLE VII. Combined predictions from supernovae signals in water and heavy water and 
corresponding models with 3 and 4 flavour mixing that are consistent with them. By "None" we 
mean none of the models of mixing that we have considered here. The notation is the same as in 
the earlier table with the ratios Rd and Rh referring to heavy water and water respectively. 







Number of events on 


deuteron 




#4/0 


sin 2 2uj 


No mixing 


3-flavours 




4-flavours 


0.960 


374 


374 




274 


0.73 


0.760 


374 


374 




281 


0.75 


0.006 


374 


374 




293 


0.78 



TABLE VIII. NC events on deuteron with E v > 3 MeV for different values of uj in the adiabatic 
case when e = 0.08. While the 3-flavour case is identical to the no-mixing case, as expected, the 
4-flavour case shows a depletion in events due to loss into the sterile channel. The ratio of the 
4-flavour to the no-mixing (or 3-flavour) case is shown in the last column. 
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Number of events on 


deuteron 


-R4/0 


sin 2 2u 


No mixing 


3-fiavours 


4-flavours 


0.960 


374 


374 


336 


0.90 


0.760 


374 


374 


338 


0.90 


0.006 


374 


374 


365 


0.98 



TABLE IX. The same as Table VIII, with e ~ so that the upper resonance is non-adiabatic. 
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FIGURES 



S m LSND 





8m ATM 






v 2 

r 5m s 


6m 2 s 



(a) 



(b) 



FIG. 1. The vacuum mass square differences in the 3 and 4 flavour schemes. In the 4-flavour 
scheme, v e and v s are predominantly mixed states of v\ and vi while and v T are that of v% and 
1/4 (2+2 scheme). The mixing between the lower and upper doublets has been chosen to be very 
small. Here S, ATM and lsnd stand for the solar, atmospheric and lsnd mass squared differences 
respectively. 




FIG. 2. Schematic drawing showing mass squares as functions of matter density in the 3 and 
4 flavour schemes. Resonances occur at two different regions of matter density, the lower one at 
5m? s . In the 4-flavour case the upper resonances consist of four close- lying resonances determined 
by Sml sm . 
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20 



60 



40 

E e (MeV) 

FIG. 3. v e d event rates as a function of the scattered electron energy E e , when the upper 
resonance is completely adiabatic. The solid line represents the no-mixing case. The dotted line is 
due to the effects of either 3- or 4-flavour mixing, which cannot be distinguished here. 
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0.6 



0.4 



5 = 5.4 10" 6 
sin 2 2cj = 0.006 



20 40 60 80 100 

E„ (MeV) 

FIG. 4. The Landau Zener jump probability Pl at the lower resonance as a function of the 
neutrino energy for the small angle MSW solution (SMA) values of the mass squared difference 
and mixing angle uj for the 3-flavour case. 
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FIG. 5. v e d event rates when the upper resonance is completely non-adiabatic. The solid lines 
represent the no-mixing case. The dotted and dashed lines are due to the effects of 3- and 4-flavour 
mixing. Results are shown for two different values of u when e = 1CT 4 . 




20 40 60 20 40 60 



E e (MeV) E e (MeV) 



FIG. 6. V e d event rates are shown as a function of positron energy, E e . The solid lines represent 
the no-mixing case. The dotted and dashed lines are due to the effects of 3- and 4-flavour mixing. 
Results for two different values of u> are shown. 
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FIG. 7. Total event rates (from combining the indivdual u e d and V e d rates shown in Figs. || 
and ^) are shown as a function of the electron/positron energy, E e , for two different values of oj, 
and for e = 0.08 so that the propagation is fully adiabatic. The dotted and dashed lines are due 
to the effects of 3- and 4-flavour mixing. Results from a 1 kton water detector (from V e p alone) 
are shown on the right, for comparison. 
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FIG. 8. Same as Fig. [7| but for e ~ so that non-adiabatic effects are included. Hence this is 
a combination of Figs. |5| and ||. 



33 



05 



03 



OS 



2.5 



1.5 



0.5 



n i i 1 1 ii i i i 1 1 ii i i i i 1 1 

E e ^5; f = 3 
e = 0.08 



e~0 



j i 1 1 1 1 in i i i 



~i — i 1 1 1 1 mi — i — i i 1 1 1 1 — i — i 1 1 1 1 





I I 

TP >> C . 

E e ^5, 




e = 0.08 











J I I 1 1 1 III I I I 



10 



- 3 10- 2 10- 1 
sin 2 2w 



10" s 10- 1 
sin 2 2cj 



FIG. 9. The double ratio Rd/h of the ratio of the total events observed through the detection 
of an electron (or positron) with E e > 5 MeV, from a future supernova explosion to that expected, 
in a heavy water and a water detector, shown as a function of the (12) mixing angle u. Solid and 
dashed lines correspond to adiabatic (e = 0.08) and non-adiabatic (e ~ 0) neutrino propagation at 
the upper resonance(s). The case for 3-flavour mixing is shown on the left and that for 4- flavours 
on the right. The vertical dotted lines indicate the currently favoured values of sin 2 2a; according 
to solar neutrino analysis. 
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FIG. 10. The same as fig. ^, with a high energy cut on the electrons, 5 < E e < 40 MeV, to 
decrease sensitivity to the high energy tail of the neutrino spectrum. 
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20 40 60 

E e (MeV) 



FIG. 11. The same as Fig. |3] for electron neutrinos in the no-mixing case (solid lines). Added 
are the dotted and dashed lines corresponding to the case when the supernova neutrino spectral 
temperatures are increased or decreased by 1 MeV in each time-bin. 




40 

E e (MeV) 

FIG. 12. The same as Fig. || for electron antineutrinos in the no-mixing case (solid lines). 
Added are the dotted and dashed lines corresponding to the case when the supernova neutrino 
spectral temperatures are increased or decreased by 1 MeV in each time-bin. 



35 





2.5 


X 

OS 


2 






II 


1.5 


X 




\ 

Q 


1 


Oh 




0.5 




1 i i i I nil | I | 



10" 3 10" 2 10" 1 10" 3 10- 2 10" 1 
sin22w sin 2 2cj 



FIG. 13. The same as fig. ^, but where the dashed curves now indicate the changes in the 
double ratio Rd/h due to variations in temperature in the supernova model by ±1 MeV. The 
statistical errors are also plotted but are too small to be distinguished. 
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